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Abstract 

We analyze flavor-changing-neutral-current (FCNC) effects in the b ^ s tran- 
sitions that are induced by family non-universal ^7(1)' gauge symmetries. After 
systematically developing the necessary formalism, we present a correlated analy- 
sis for the AB = 1,2 processes. We adopt a model-independent approach in which 
we only require family-universal charges for the first and second generations and 
small fermion mixing angles. We analyze the constraints on the resulting param- 
eter space from Bg — Bg mixing and the time-dependent CP asymmetries of the 
penguin-dominated B^ — > {iTjCp,!]' , p,io, fo)Ks decays. Our results indicate that 
the currently observed discrepancies in some of these modes with respect to the 
Standard Model predictions can be consistently accommodated within this general 
class of models. 
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1 Introduction 



The origin of CP violation, which was first observed in the kaon system four decades ago 
[T], has remained one of the fundamental questions of elementary particle physics. In 
recent years, the B factories have established that the Standard Model (SM) picture of 
CP violation, in which all CP-violating effects are generated by the single phase (5ckm 
in the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix [21 |3], is consistent 
with the observed pattern of CP-violating phenomena in both the and K meson 
systems [1]. However, as the SM cannot account for the baryon asymmetry in the 
Universe todayj^, new physics (NP) is necessarily required to describe all observed 
phenomena with CP-violation involved. 

One arena to seek the NP effects is in flavor-changing neutral current transitions 
(FCNC) where the SM contributions flrst appear at the one-loop level and the NP ef- 
fects can be competitive. The emblematic set of such processes is the set of 6 — > s 
transitions, which include Bg — Bg mixing and the set of neutral Bd meson decays which 
occur via b sqq (g = u, d, c, s) transitions. Several of these processes are also of 
interest because recent measurements exhibit discrepancies with the SM predictions at 
the level of a few standard deviations, which may suggest the intriguing possibility of 
physics beyond the SM. The current status of the data is as follows: 

• Bg — Bg mixing phase. The standard way to parametrize NP in Bg — Bg mixing is 
to express the off-diagonal mixing matrix element as follows: 

= (Mit)sMCB.e2^<^Br. (1.1) 

The SM predicts that C^^ = 1 and = 0. Though the data indicate that Cb^ does 
not differ signiflcantly from unity, the results of a recent analysis [0] suggest that 4>^^ 
deviates from zero at the Sex level (see Table [T]). This analysis combines all the available 
experimental results on Bg mixing, including the new tagged analyses of Bg ipcj) by 
CDF [7] and D0 [8J (note that no single measurement yet has a Sa signiflcance.). The 
discrepancy disfavors NP scenarios which obey minimal flavor violation (MFV), i.e., 
with (f)^ ^ 0, and instead suggests NP which exhibits flavor violation in the h ^ s 
transitions (e.g., see [9] and references therein). For convenience, in Table [l] we also give 
the data in terms of A^^ /A^J^ and which are related to Cb^ and according to 

C,.e-^^.r = 1 + 4^e-^?^ (1.2) 



• CP asymmetries in neutral Bd decays. The set of neutral Bd decays in question 
is the set of QCD penguin-dominated charmless decays that occur via h sqq (g = 
ti, (i, c, s) transitions. The CP asymmetries of such decays into a flnal CP-eigenstate 
fcp are given by 

V{Bd{t) ^ fcp) - T{Bd{t) ^ fcp) 



T{Bd{t) ^ fcp) + T{Bd{t) ^ fcp) 
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Observable la C.L. 2a C.L. 



<n (SI) 


-20.3 ± 5.3 


[-30. 


5,-9.9] 


(S2) 


-68.0 ± 4.8 


[-77.1 


B,-58.2] 




1.00 ± 0.20 


[0.6i 


B,1.51] 



0fP[°] (SI) -56.3 ± 8.3 [-69.8,-36.0] 

A^P/ASM (SI) 0.66 ± 0.28 [0.24,1.11] 

0fP[°] (S2) -79.1 ± 2.6 [-84.0,-72.8] 

A^^/A^^ (S2) 1.78 ± 0.03 [1.53,2.19] 



Table 1: Fit results for the Bg — Bs mixing parameters [6 . The two (j)^ solutions 
("SI" and "S2") result from measurement ambiguities; see |6j for details. 

= -Cf^, cos(AMB,t) + Sf^, sin{/\MB,t), (1.3) 

in which Cf^p and Sf^p are direct and mixing-induced CP asymmetry parameters. 
The SM predictions for many decays of this type, including ipKs and Bd — > 

(0, 1]', Tc, p, uj, fo)Ks are as follows: 

- VfcpSfcP = sin 2/5 + 0{\'), Cf,p = + ^(A^), (1.4) 

with (3 = aJ:g[—{VcdV*fJ)/{VtdV^l)], X = sinOc being the Cabibbo angle, and rjf^p = 
±1 being the CP eigenvalue for the final state fcp- However, the central values of 
sin 2/3 directly measured from the penguin-dominated modes are systematically below 
the SM prediction and the results obtained from measuring the charmed B^ ipKs 
mode. Meanwhile, the central values of the direct CP asymmetry measured from B^ 
(f)Ks and B^, — >■ uKs modes are also small compared to that obtained from the B^, 
ipKs mode (see Table [2]). Given that the B^ decay is dominated by tree-level 

amplitude in the SM, large absolute values for ASf^p = —Vfcp^fcp + VtpKs^ipKs ^-^d 
^^fcp ~ ^fcp ~ ^ipKs ni^y imply interesting NP in the b s transitions. . 

To account for these discrepancies appearing in Bg — Bg mixing and B^ decays, a 
number of NP scenarios have been studied, including low energy supersymmetry and 
models with warped extra dimensions, among others [TT]. In many of these scenarios, 
the effects of NP in the h ^ s transitions are loop-suppressed and can compete with SM 
contributions. The most popular and well-studied scenarios are models with minimal 
flavor violation (MFV), in which the only source of CP violation is the single irremovable 
phase of the Cabibbo-Kobayashi-Maskawa mixing matrix. MFV scenarios, however, face 
difficulties in that they do not generally allow for a nonvanishing 0^^. 
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fcp 


-VcpSfcp (If^ C.L.) 


C;^,(laC.L.) 


ijKs 


+0.672 ± 0.024 


+0.005 + 0.019 


4>Ks 


, Q 44+0.17 


-0.23 + 0.15 


v'Ks 


+0.59 ±0.07 


-0.05 + 0.05 


TiKs 


+0.57 + 0.17 


+0.01 + 0.10 


pKs 


+0.63l°il 


-0.01 + 0.20 


LoKs 


+0.45 + 0.24 


-0.32 + 0.17 


foKs 


+0.62l°i^ 


0.10 + 0.13 



Table 2: World averages of the experimental results for the CP asymmetries in Bd 
decays via h qqs transitions [10] . 

In this work, we will study the constraints from b s transitions on models with 
family non-universal gauged f/(l)' symmetries. Additional f/(l)' gauge symmetries are 
present in many well-motivated extensions of the SM, such as grand unified and/or string 
models (e.g., see [12\ for a review). Such scenarios are of particular interest because 
unlike the scenarios studied above, they allow for the intriguing possibility of tree- 
level FCNC, with contributions that are competitive with the SM even for small f/(l)' 
couplings. Depending on the details of the model, family-dependent U{iy scenarios 
can result in new FC operators and/or modified Wilson coefficients to the existing SM 
operators in the operator product expansion, providing a rich framework beyond MFV 
to explore FCNC and CP-violating effects. 

We follow the general framework for addressing Z'-induced FCNC as developed 
in [13] and systemize its application to b s transitions. Rather than considering 
specific U{iy models, we adopt a model-independent approach in which the main re- 
strictions are family universal charges for the first and second generations and small 
fermion mixing angles. We also neglect the effects of Z — Z' mixing (which are known to 
be small), and assume the absence of any exotic fermions that could mix with the usual 
SM fermions through non- universal Z' couplings, which may also result in nontrivial 
FCNC effects (e.g., see [12]). 

This work is an extension of our earlier work [T^ [TB] , in which we performed a cor- 
related analysis of the AB = 1,2 processes mentioned above for a specific set of U{iy 
scenarios. That analysis was in contrast to other studies of f/(l)' scenarios based on 
mode-by-mode analyses [15]. The purpose of this paper is twofold: first, to provide 
more details of the formalism and analysis than were given explicitly in [H], and sec- 
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ond, to analyze a more general set of f/(l)' models. Our results demonstrate that the 
6 — > s transitions not only place important constraints on family non-universal Z' cou- 
plings and mass scale, but also that family non-universal f/(l)' scenarios can explain the 
currently observed discrepancies with the SM predictions for — Bg mixing and the 
time- dependent CP asymmetries of the penguin-dominated Bd {n, (f),ri' , p,uj, fo)Ks 
decays. 

This paper is structured as follows. We begin by providing an overview of the 
formalism of the Z' induced FCNC effects in the b —>■ s transitions and present the 
effective Hamiltonian for the processes of interest at the b quark mass scale in Section 2. 
In Section 3, first we analyze the FCNC constraints within several special limits of the 
general f/(l)' parameter space, and then turn to a more general analysis. Our summary 
and conclusions are presented in Section 4. 



2 Theoretical Background 

2.1 Formalism of Z'-induced FCNC Effects 

The general framework for studying Z'- induced FCNC Effects has been developed in |13j . 
In this section, we will systematically formalize its applications to the case of the b ^ s 
transitions (the generalization to b —>■ d transitions is straightforward). 

We begin by considering the SM extended by a single additional f/(l)' gauge sym- 
metry (the generalization to multiple U{1)' gauge symmetries is straightforward). In 
this theory, the neutral current Lagrangian in the SM gauge eigenstate basis is given by 

Cnc = -eJLK - giJzZ, - g2rz'Z',, (2.5) 

in which is the f/(l)„,, gauge boson, is the massive electroweak (EW) neutral gauge 
boson, Z'^ is the gauge boson associated with the additional Abelian gauge symmetry, 
and gi = g/ cos 9w and g2 are the gauge couplings of the Z^ and Z'^ bosons, respectively. 
The currents are given by 

= EE^^^'^^^^+^^^«]^- (2-6) 

■(/> i 

Jz' = + (2-7) 

in which ip labels the SM fermions, i and j are family indices, and Pr^l = (1 =t 75)/2- 
The (family universal) SM chiral charges are given by 

ef^ = tt^ - sin^ ewQi.,, ef^ = - sin^ dwQ^^, (2.8) 

in which t^^ denotes the third component of the weak isospin and are the electric 

charges of ipL,R- Without loss of generality, the Z' chiral charges can be diagonalized by 
choosing the appropriate gauge basis for the fermions: 

itr = (2.9) 
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In particular, SU (2) l symmetry requires that 

e7-^et, eT^eT- (2.10) 

If the diagonal U{iy chiral charges are non-universal, flavor- changing (FC) Z' couplings 
are generically induced by fermion mixing. The fermion Yukawa matrices in the 
weak eigenstate basis are diagonalized by the unitary matrices V^^ such that 

and the CKM matrix is given by 

V^CKM = Vu^Vl. (2.12) 
Hence, the chiral Z' couplings in the fermion mass eigenstate basis take the form: 

B^- = V^,^-Vl , B^- ^ V^J^-Vl . (2.13) 

However, it is known that the constraints from K — K mixing and from — e conversion 
in muonic atoms exclude significant non-universal effects for the first two families, which 
suggests that 

/ Bti''' \ 

= Sfi^'^ Bti'" , (2.14) 

\ T3i'L,R* ryi'L.R* TDi>L,R j 

\-Dl3 -D23 -^33 / 

at least for the down-type quarks and e, fj,, r leptons. The most straightforward way to 
achieve this coupling structure is to assume universal 11(1)' charges for the down-type 
fermions of the first two families, i.e., 

nt''" 

e'^^-^ = e?^'« I . (2.15) 

V e; 

With the unitary matrices V^^ ^ written as 

^^.,H= f y''^'" (2-16) 
where W^^ ^ is a 2 x 2 submatrix, one obtains 



B^ 



(2.17) 



6 



Therefore, in the hmit of small fermion mixing angles or small X^^ ^, Y^^ ^ elements, a 
Z' coupling structure of the type given in Eq. (2.14) is produced, in which 



^11 = ^1 , 

^13 5 ^23 ~ 



-Dq, — e 



'33 

0(X 



~^L,R 

3 



(2.18) 



such that 5^3^'^ and st^''' are in general both complex parameters. 

EW symmetry breaking induces Z — Z' mixing, such that the gauge eigenstates Z^ 
and Z'^ are related to the mass eigenstates zj^^^ (n = 1, 2) by an orthogonal transforma- 
tion. In the mass eigenstate basis, the Lagrangian couplings are given by 



*-NC 



[gi cosOJ^ + 92 sinej^,] Z^'^ - [-g, sinOJ^ + g2 cosOJ^,] Zf, (2.19) 



where Q is the Z — Z' mixing angle, is given in Eq. (2.6), and J^, is of the form 
of Eq. (2.7) with e^^^R replaced by B'^^'R from Eq. (2.13). In this analysis, we neglect 
kinetic mixing since it simply amounts to a redefinition of the unknown Z' couplingsj^ 
At the EW scale, the tree-level four-fermion interactions are described by the product 
of gauge currents 



{PeffJz'^ + 2wJz ■ Jz' + yJz''^) 



-4Gp 



V2 



\^ mn'-^mn ~ ^mn'^mn ~ mn mn ~ mn mn 



. (2.20) 



In Eq. (2.20), the local current-current operators ar^(i,j, m, n are family indices): 

= {i^^rPLi^i) iXml.PRXn) , f^n = {i^il'PRi^,) {Xml.PLXn) , (2.21) 

and the coefficients are 

C/^ = n ffS S e'^^e^^ + w8 -e^^B^^ + e^^ B"^^ + vB"^^ B^^ 

C/^ = n ffS S e'^^e^^ + w^ -e'^^B^^ + e^^B"^^ + vB'^^B^'^ 
o ff?i- ?i e'^^e^^ + w?i -e'^^B^^ B"^^ -\- vB"^^ B^^ 

n ffS S e'^^e^^ + w^ -e'^^B^^ + e^^ B"^^ + vB"^^ B^^ 



mn 



mn 



(2.22) 



^Kinetic mixing allows the redefined Z' charges to have a component of weak hypercharge, which 
would otherwise not be allowed. This feature is irrelevant for the purposes of this paper. 

^These operators are not all independent. For couplings of four fermions of the same type, ip = Xj 
e.g. four charged leptons, one has Sli„ = 5^^", S^^ = S^"" and T^„ = 
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in which 



Peff 



W 



Pi cos^ + p2 sin^ 9, 



pa 



Ml cos2 e 



w 



92 



9i 



sm9 cos9{pi — P2) 



(pi sin^ 9 + P2 cos^ 



(2.23) 



In Eqs. (2.23), Ma denotes the masses of the neutral gauge boson mass eigenstates, 
and 9w is the EW mixing an gle. W e do not specify the ip and x dependence of the 
coefficients C, C, D, D in Eqs. (2.22), which can be understood from the context. 

For b ^ s transitions, the local operators are given by S'^„, 5'^„, T^^^ and T^^, with 
coefficients that are given by 



bs 
mn 



/-Ids 

'pit 

r 



XL 
mn 1 



bs 
mn 



D 
D 



bs 



bs 



WOmntm ^bs ^ iJ^bs^mn ) 

wd^ne'^^^Bg + yB^^B 



mn 



(2.24) 



With the Z — Z' mixing angle neglected, the coefficients can be written as 



C: 



bs 
mn 



C: 



bs 
mn 



D 



bs 



D 



bs 



f g2Mz 

\giMz' 

( g2Mz 
\giMz' 

( g2Mz 
\giMz' 



^bs-^mn 1 



-^bs-^mn ' 



nL T^XR 

bs mn 1 



uR nXL 



(2.25) 



For convenience, in the following we will resolve the factor g2Mz/{giMz>) into the B 
elements or the chiral couplings. 

At tree level, there are three classes of 6 — s transitions which are sensitive to the 
possible NP effects that result from an additional family non-universal U{iy symmetry: 
b sqq transitions, b — >■ sll transitions, and Bg — Bg mixing. Here "q" and "1" denote 
quarks and leptons, respectively. For the b sqq transitions, the Z' effects are described 
by the effective Hamiltonian 



H^^ib ^ sqq) 



2Gf 



'sb)v-A J2K,im)v-A + Dll{qq)v+A) 
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+ {sb)v+A Y.^D':My-^ + Climv+A)) + h.c, (2.26) 



in which the sum is over the active quarks for a given process. These Z'-induced FCNC 
effects can be understood as corrections to the SM operators or to the new penguin 
operators defined in [Appendix A[ since both lead to the same hadronic matrix elements. 



Explicitly, comparing Eq. (2.26) with 

sqq) = 



n/Z 



V2 

+ {sh)v+A 



V+A 



~ ~ 3 ~ ~ 3 

( (AGs + AC^-eq){qq)v+A + {AC, + ACj-eq){qq)v-A 



+ h.c., 
(2.27) 



results in Auq equations {uq is the number of active quarks in the final states): 





+ AC9-e, = 


T/, T/* QQ 

^tb Vts 


AC, 


+ ACr^eg = 


^tb yts 


ACs 


~ 3 

+ ACg-eq = 


2 ^bs 
Vtb Vts 


AC, 


+ ACr-eq = 


-2 ^bs 

Vtb Vts 



(2.28) 



where AC denotes Z' correction to the Wilson coefficients of the SM operators and AC 
denotes the Wilson coefficients of the operators beyond the SM ones. For charmless 
processes with q from the first two families, these equations are solvable because of the 
following relation obeyed by the down-type quark couphngs: 



B 
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B. 



22 1 



(2.29) 



which is extracted from Eq. (2.14). The Z' corrections to the Wilson coefficients are 
then found to b43 



AC, 



{Ct + 2C^^) 



AC 
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tbVts 



( /^bs /^bs\ 



■^Though the solutions to Eq. (2.28 1 are not unique in the case with = 1, the physics is unaffected 
since it is only sensitive to the linear combinations on the left hand side of Eq. (2.28). For the charmed 



processes where generally we have Ug — 1 , the formula in Eq. (2.301 can also be applied as long as 
Eq. (2.291 holds for the up-type quarks. If Eq. (2.291 does not 



lold, then the "uu" indices in these 



formula need to be replaced by "cc" . 
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AC, 
AC,-- 



2 
2 



7 [Dt + '^D%) 



ACj 
AC, 
ACj 



We pause here to comment on subtleties in Eq. (2.30). Recall that in the limit of 



4 


\ uu 


- D^d) , 








4 


y-'uu 






mbV^s 




4 


( 

\ uu 




(2.30) 





small fermion mixing angles, Eq. (2.18) holds for the down-type quarks. To obtain 



the CKM matrix as given in Eq. (2.12) without requiring fine-tuned cancellations, the 



mixing angles for the up-type left-chiral quarks should also be small in this limit. Due 



to the SU{2)l constraint of Eq. (2.10), therefore, Eq. (2.18) can also be applied to the 
up-type left-chiral quarks. In this case, it is straightforward to see that 



Bi 



B. 



dd 



and hence 



AC, 



-c: 



bs 



- ei = 0, 



ACr 
AC, 



(2.31) 



tbVts 



VtbVts 



-D 



bs 
dd- 



(2.32) 



Note that a relation similar to Eq.(2.10) does not exist for the right-chiral SM fermions, 
so ACj and ACg are generically non-trivial. In regards to the color-allowed penguin 
operators, their Wilson coefficients are corrected by Z' effects only at the loop level 
where the color-indices are mixed by gluons. Since these effects suffer loop and Z' mass 
double suppressions, we will not consider them further in this paper. 

For the b — > sll transitions, the Z' contributions to the effective Hamiltonian are 



n^sib - sll) 



2Gf 



7sh)v-A{Cli\ll)v^A + Dlt{ll)v+A 



Comparing Eq. (2.33) with 



HffAh^sll) 



+{sh)v+A{Dl^{ll)v-A + Cli\ll)v+A + h.c. 



Cp 



(2.33) 



'^tb^ts ( ACijvQw + ACiqaQwa 



+ACwQw + ACioaQwa + h.c. 



(2.34) 



one can see that the Z' corrections to the Wilson coefficients take the following form: 

2 



ACioA 



VtbV,s 
2 



VtbV,l 



-{ct+D^n, 
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Vtb Vts 

^C^oA = -T^M'~D',n. (2.35) 

Vtb Vts 



Note that if the leptons in the process are neutrinos, Eqs. (2.35) reduces to 

VtbVts 
VtbVts 
VtbVts 

^C,oA = T^Mu (2.36) 

since right-handed neutrinos are generally decoupled at low energy scales. 

For Bg — Bs mixing, the Z' corrections to the effective Hamiltonian take the form 

n!;,{Bs-Bs) = ^(c,^:(.-6)v-A(s6)y-A + /^L^(^6)y_A(s-&)y+A 

+D',l{sb)v+A{sb)v-A + C',:{sb)v+A{sb)v+A) + h.c. (2.37) 
Once again, upon comparing this expression to 

n^ffiBs -Bs) = -^ (ACf ^gf^ + ACf + 2ACi^Qi^^ + h.c, (2.38) 

the Z' corrections to the Wilson coefficients are easily determined to be 

ACf= = -Cts, 

Acf^ = -ct: 

AC^^ = -l[Dt + &:;) = -bt. (2.39) 

As in the b —>■ sqq transitions, the Z' effects only correct the Wilson coefficients of the 
color-allowed operators at a higher loop level, so we will not consider them further here. 

To summarize, in Table ([3]) we classify the tree-level Z' contributions to the b ^ s 
transitions according to whether they are relevant or irrelevant to the SM operators. 

Before considering the general parameter space, it is worthwhile to consider a few 
special limits: (1) the LR limit: = = (2) the LL limit: e'^^ oc J; and 

(3) the RR limit: e^^ oc /, where / is the identity. The Z' corrections to the Wilson 
coefficients in these limits are summarized as follows: 



11 



Table 3: Classification of the tree-level Z' corrections to the Wilson coefficients in 
the 6 — > s transitions. 





SM operators 


Beyond SM Operators 


h — > sqq 


AC„, a = 3,5,7 


ACa, a = 3, 5, 9 


b^sll 


ACa, a = 9V, lOA 


ACa, a = W, lOA 


Bg — Bg mixing 




AC'f^a = l,3 



(1) LR limit: B^^ = B^ 



bs- 





ACf ^ = 




ACs = 


AC,^- 


T/ ^r^^bs^ddi 

nb^ts 


ACa = 


AC5 = - 


L R R 


ACy = 


ACg = - 


TjL ( jjR jjR\ 

'^nbVts 


ACgv = 


ACgy = 


2 L L 

..*Bf,g{Bii +Bii), 
nbVts 


ACloA — 


ACloA = 


2 

nbVts 



(2) LL limit: e'^« oc / 



(3) RR limit: e'^^ oc I 



ACf^ 
ACa 

AC, 

AC-! 

ACgv 

ACioA 

AC-f^ 



2 



pL pL 
^bs^ddi 



VtbV^s 
2 

4 

2 L L 

nbVts 
2 



5,^, + 25,^,) , 

pL /pR p_RN 
^bs [Buu - Bdd) , 



VuMtg 



-Bt{-B^i-rBi^). 



(2.40) 



(2.41) 
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2 
2 

VtbVts 
^tbVts 

ACioA = -TT^B^,{-B[i + Bll). (2.42) 

We will focus on the correlations between Bg — Bs mixing and the hadronic B^ meson 
decays. For the latter, though Z'-mediated effects can occur in both the QCD and EW 
penguins, we make a conservative assumption in this paper that they are mainly manifest 
in the EW penguins, such that |AC3,5| -C [ACyj, as suggested in [TTl [15]. With this 
restriction, there are only three relevant parameters for each special limit: the modulus 
of 5^ (or -B^), its phase 0^ (or 0^), and the real -Bj^l— ~-^mL/2)- These parameters 
need to satisfy 



BL\ < \B!i,\ « \Bf,\ (2.43) 



in the LR and LL limits, and 



\B^s\ < \Bfdl \B^d\ « \B!,\ (2.44) 

in the RR limit. Here < l-B^j^l is due to the fact that, under the assumption 

of small fermion mixing angles, the modulous of off-diagonal elements in the coupling 



matrix in Eq. (2.14) should be smaller than that of diagonal ones; l-B^^I <^ is due 



to lACs^sl <^ IAC7I. Later in the paper, we will consider the more general parameter 
space for the Z'-mediated effects in the EW penguins, which has five free parameters: 
\B;:;% 4f and Bf,. 

2.2 Effective Couplings at the b Mass Scale 

To achieve sufficient precision for these observables, it is necessary to have an accurate 
knowledge of the relevant Wilson coefficients at the b quark mass scale rrib = 4.2 GeV. 
The Wilson coefficients at the b mass scale can be obtained as follows: 

C{m,) = U{m,, Mw) C{Mw), (2.45) 

where C is a vector with entries consisting of the Wilson coefficients and U is the evolu- 
tion matrix. The observables can then be expressed in terms of the Wilson coefficients 
at the rrib scale (for general discussions, see e.g. [H]). All parameter values used in our 



calculations are summarized in Appendix C 
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• Bs mixing. Following [18], the NP probes Cb^ and 0^^, which are defined in Eq. ( 1.1 ), 
are calculated to be 

Cb^c^'^Bs = I _ 3,59 X 10^(ACf= + ACf^ + 2.04 x 10^ AC^^ (2.46) 

at the nib scale. The large coefficients of the correction terms are due to the fact that 
the NP is introduced at tree- level while the SM limit is a loop-level effect. 

• Bd ^ ttKs decays. The Bd t^Ks decays have recently received considerable 
interest in the literature (see e.g. [HI [191 EDI EH ESI [23] ) . In [H], it is pointed out that 
a deviation of St,Ks from its SM value can be understood as a modification of the ratio 

in which T, C and P denote the color-allowed tree, color-suppressed tree, and EW pen- 
guin contributions in the decay amplitude, respectively. In the class of models considered 
here, the family non-universal Z' interactions modify ge**^ through the relation 

qe^'t' = 0.76(1 + 158. 1AC7-102.4AC'9), (2.48) 

in which q and are given by 0.76 and zero, respectively, in the SM limit. 

• Bd ^ {ip,(j),r]', p,uj, f^)Ks decays. The direct and the mixing-induced CP asymme- 
tries in the Bd hadronic decays are parametrized as follows: 

1- l^/cpP c _ 2Im[A/^/ 



In the above, A/^^p is defined by 



A/cP--^/c.^^, (2.50) 

PBd ^fcp 



with 



Ars,=0 |(MbJi2| 

Here qs^ and are Bd mixing coefficients 

\Bl) = pB^Bd) + qB^Bd) 

\Bh) = pB,\Bd)-qB,\Bd), (2.52) 

Mb^ is the Bd — Bd mass matrix, and Af^p is the decay amplitude of Bd fcp (^/cp 
is its CP conjugate.). 
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The SM predicts that = P = arg [—(ycdV*f^) / {VtdV^l)] and that a non-trivial weak 
phase enters Af^p only at order C(A^). Therefore, for time-dependent decays proceeding 
via b — > sqq{q = u,d,c,s), including Ba — > ipKs and penguin-dominated modes such 



as Bfi — > {(f),-)]' ,71, p, uj,f^)Ks, the relations in Eq. (1.4) are obtained. However, these 
results are greatly changed with the involvement of family non-universal Z' bosons, since 
this allows for a new weak phase to enter Aj^p at tree level. Following Ali et. al. |24j . 
the Xfcp parameters of Ed —>■ {ip, 4>, r]', vr, p, u, P)Ks are given by 

A^;^^ = (-0.63 + 0.74i) (2.53) 
1 + (0.18 - 0.01z)(AC7 + AC-jY - (0.06 - 0.04z)(AC9 + ACg)* 
1 + (0.17 + 0.01i)(AC7 + ACr) - (0.05 + 0.05i)(AC9 + ACg) ' 

A^^^ = (-0.70 + 0.70i) (2.54) 
1 + (14.57 + 5.88z)(AC7 + ACy)* + (15.08 + 5.92i)(AC9 + ACg)* 
1 + (14.39 + 5.64i)(AC7 + ACy) + (14.90 + 5.66i)(AC9 + ACg) ' 

A^,;^^ = (-0.70 + 0.69?) (2.55) 
1 + (2.11 + 0.67i)(AC7 + ACr)* + (2.10 + 0.54i)(AC9 + ACg)* 
1 + (2.08 + 0.65i) {/\Cj + ACr) + (2.07 + 0.52i) ( ACg + ACg) ' 

= (-0.74 + 0.66i) (2.56) 

9) ' 



1 - (38.75 + 3.29i)(AC7 + ACr)* - (47.95 + 4.1l2)(AC9 + AG 
1 - (38.11 + 5.23i){AC7 + ACy) - (47.15 + 6.5Qi){AC9 + AG 

(-0.71 + 0.70i) ( 
1 + (31.97 + 4.76i)(AC7 + AC7)* + (18.84 + 2.75z)(AC9 + ACg)* 
1 + (31.81 + 4.67i)(AC7 + AC7) + (18.74 + 2.70z)(AC9 + ACg) ' 



XuKs = (-0.71 + 0.70i) (2.57) 

' 

A^o^^ = (-0.70 + 0.702) (2.58) 
1 + (3.19 + 0.93i)(AC7 + AC7)* - (0.12 + 0.15z)(AC9 + ACg)* 
1 + (3.16 + 0.90i)(AC7 + AC'7) - (0.12 + 0.15i)(AC9 + ACg) 

In contrast to the — > ipKs decay, in which the NP effects are suppressed by the 
SM tree-level contribution, family non-universal U{1)' couplings indeed result in sizable 
corrections to A/^p for the penguin-dominated modes. 
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Figure 1: Correlated constraints on iBf^J and are presented. In these two panels, 
random values for and (/>§^ from the experimentally allowed regions 
(see Table 1 1 are mapped to the \B^J 



£, plane using Eq. (3.59), with an 



assumed 25% uncertainty (a typical value from non-perturbative effects) 
assumed for the coefficients. The left (right) panel is the LR (LL) limit. 



3 Results and Analysis 

3.1 Correlated Analysis (I) — Special Limits 

In this section, we will present a correlated analysis of the AB = 1,2 processes which 
occur via b ^ s transitions, focusing first on the special limits of the parameter space 



as presented in Section As the physics of the RR limit is very similar to that of the 
LL limit, we focus in this paper on the LR and LL limits as representative examples. 

We first consider the constraints on this class of family non-universal U{iy scenarios 
which arise from Bg — Bg mixing. With the renormalization scale chosen as the 6-quark 
mass, rrih = 4.2 GeV, the NP probes Cb^ and are given by 

C^^e^^'^^r = 1 + 1.32 X lO^ACf^ 

Cs^e^^^^r = 1 _ 3.59 X lO^ACf (3.59) 



in the LR and LL limits, respectively. These conditions involve two of the three free 

^bs\ 



parameters of each limit: and The experimental constraints on these param- 



eters from Bg — Bg mixing are illustrated in Fig. [T] The left panel corresponds to the 
LR limit and the right one corresponds to the LL limit; in this section we will present 
the results for these two limits together, so that it is easy to make comparisons between 
the two cases. In each case, there are two separate shaded regions, corresponding to the 
two solutions (see Table [T). For each region, the various colors of the points specify 
the different confidence levels (C.L.) of the relevant Cb^ and values. To explain the 
observed discrepancy in Bg — Bg mixing from the SM prediction, is required to be 
~ 10~^. This reflects two facts: (1) unlike (ps^, the modulus Cb^ does not deviate from 
its SM prediction signiflcantly (the experimental value of Cb^ has an at most 0{1) shift 
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Figure 2: The constraints on from g'e*'^ are sliown. Tlie points from the {Bj^J—cp^g 
plane (see Fig. [l]) are randomly combined with scattered points of 
(10~^ < < 10~^) and then mapped to the gcost/) — qsintfi plane 

according to Eq. ( 3.64[ ). The colors of the points in this plane indicate 
the C.L. that their inverse images represent in Fig. [T] The two dashed 
lines specify the experimentally allowed ranges that result from the 
fit of the B ^ ttK (and B vrvr) data at la and 90%(~ 1.7a) C.L., 
respectively |19| . The left (right) panels are the LR (LL) limits. 



from its SM prediction at 2a C.L.); (2) the Z' corrections are from tree level, and hence 
can easily explain this small deviation (only a small coupling is necessary, according to 
Eq. (3.59) and Eq. (2.40), Eq. (2.41)). The smallness of is generically consistent 



with our assumption of small fermion mixing angles, since is proportional to them 



as well as to g2Mz/{giMz') (see Eq. (2.18) and the comments under Eq. (2.25)). In 
addition, due to the smallness of l-B^I, the experimental constraints from the branching 
ratio By{Bs — > /x"*"//") can be easily satisfied. For details, see Appendix B 

Before move to the other 6 — s processes, we have some comments on the influence 
of Bg — Bg mixing on another FCNC Z' coupling B^f. In the SM, the mass differences 
of B^ and Bg masons are predicted to be (e.g., see [25] ) 



AM: 



SM 



(0.53 ±0.02) 



(19.3 ±0.6) 



td\ 



Br, 



0.0082/ V200MeVy 0.85 

2 / r \ 2 



B 



ps" 



ts\ 



0.00405 J V240MeV J 0.85 



B 



:PS 



(3.60) 



with fBd,Bs being decay constant of B^^Bs) and B being bag factor. Comparing with 
the experimental data [26j 



AM, 
AM. 



AM^P 

^Ml^il + ^^) = 0.507 ± 0.005ps-^ 



AMf^(l ± 



AM; 



NP 



AMSM' 



17.77 ±0.12ps-\ 



(3.61) 
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Figure 3: The time-dependent CP asymmetries of the charmed — > ij^Kg decay are 
presented (with \Vub\ = 3.51 x 10~^, as is used in the SM calculation |28]). 
The box is at la C.L. and the dark point is the SM limit. The left (right) 
panels are the LR (LL) limits. 



we have the relation 



(3.62) 



with 



AMP 



\B 



L,R\2 
bd I 



AM|^ 
AMSM 



^ 0.04. 



(3.63) 



of -B^'^, a fine-tuning of 0(10%) level therefore is needed in l-B^'"] to satisfy the ex- 
perimental constraints from AM^. In this paper, we will work under the assumption of 
negligible |, and therefore will neglect its possible effect in the NP observables. 

The second process of interest is tcKs- The time-dependent CP asymmetries 

of these decays can be sizably affected by NP, as has been pointed out in [T7]. The 
experimental constraints on qe^'^ (defined in Eq. (2.47)) for different C.L.'s from the 
B — >■ ttK (and the B —>■ vrvr) data have previously been obtained in [H]. In Fig. |2| we 
illustrate how is constrained through ge**^, using the parameter values of and 



But, according to Eq. (|2.18l), the modulous of Bh^ usually is comparable with that 



obtained in Fig. [T| along with the following relations: 



qe 
qe 



0.76(1 + 55.7AC7) 
0.76(1 + I58.IAC7) 



(3.64) 



which are valid in the LR and LL limits, respectively. There are two distribution regions 
which are specified by different colors in each panel, again due to the two solutions. 
Note that in the LL limit, the shaded region passes through both of the minimal points 
that were found in the fit of the B — > ttK and B ^ tttt data in 
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In the scenarios under consideration, the constraints from Bg — Bg mixing and Bd — >■ 
TcKs decays place bounds on each of the three free parameters |-B^|, and B^i- The 
natural question is then whether the experimentally allowed values for these parameters 
also satisfy the constraints resulting from the possibly anomalous values of AC/^p and 
ASf^p in the remaining penguin-dominated Bd —>■ (0, f]', p, uj, /o)-R's decays. To address 
this issue, we assume a 15% uncertainty in the SM calculations for each of these modes 
(as well as for the Bd t^Ks mode) and a 25% uncertainty for the NP contributions. 
Here 15% is a typical uncertainty level for the hadronic matrix elements of the SM FC 
operators (e.g., see [27]) and is also the least necessary one to explain the experimental 
data of C^Ks cind S^Ks iii the SM (see Fig. [s] where the NP effects are negligible). As 
for the difference of the uncertainty levels between the SM and NP calculations, it is 
caused by the fact that the hadronic matrix elements of the FC operators in the SM 
are better understood than they are for the NP operators. In Fig. [4], we systematically 
illustrate the time-dependent CP asymmetries in the penguin-dominated modes, using 



the parameter values obtained above and Eqs. (2.53)-(2.58). For these modes with the 
exception of Bd — >■ pKs, there are 0.5 ~ 2a deviations for C/^p, Sf^p, or both. Though 
our model only induces negligible effects on Bd — Bd mixing under the assumption of 
small l-Bftdl, due to the interference effects betwe en th e Bd — Bd mixing phase and 0^ 

which affects the decay asymmetries in Eq- (2.50), the points in Fig.|4jare scattered 

away from the SM limits. This results in a dispersion such that there are always some 
points lying in the la region for each of these modes. 

To show that all of the constraints can be satisfied simultaneously, we have carried 
out a correlated analysis among the Bg — Bg mixing and the Bd — > (tt, 0, p, a;, fo)Ks 
CP asymmetries. The allowed values for |-B^|, and B^^ in the LR and LL limits 
are illustrated in Fig. [5j We see that in Fig. |5] there indeed exist parameter regions 
where the anomalies in Bg — Bg mixing and the time-dependent CP asymmetries of 
Bd —>■ (tt, 0, rj', p, u, fo)Ks, can be explained by NP at reasonable C.L.. The allowed 
and 0^ values can explain both solutions of Bg — Bg mixing phase; and the allowed |-Bj^| 
values vary from 0.08 to smaller values. If we want to get a better fit for C(^-„^if,^ri',p,ojjo)Ks 
and 5(7r,</,,r?',p,w,/o)-fs' \Bdd\ ~ 10^^ is typically required. To see this point, we take for 
example Ct^Ks and S^Ks ; and map the points in Fig. [5] back to the q cos — g sin plane, 
as illustrated in Figjo] In this figure we see that the points with B^^ < —0.01 are closer 
to the minima of the fit of the B — >■ ttK and B ^ hit data, leading to a better fit 
compared to the one obtained in the SM limit. Therefore, |-Bj^| ^ 10~^ is important in 
improving the agreement with experimental data in the penguin-dominated Bd decays]^ 

The favored parameter values for |-B^| are interesting for collider detection. For 
iydpe'^''Vd^)ii ~ C»(l), this implies that giMz>/{g2Mz) ~ 10 - 100 or a TeV scale Z' 
boson for g2 ^ gi, a range approachable at the LHC (e.g., see [HI [29]). This fact is also 



important for the effective Lagrangian in Eq. (2.20) which is obtained by integrating 



out the Z' boson. While applying it to our analysis, we neglected the effects of the 



'*This effect can also be seen by requiring a smaller C.L. for the fit of C(^c/>,7i' ,p,ujjo)Ks ^^'^ 
S((p,r]' ,p,uj,fo)Ksi which has been shown in the LR limit in Fig. 4 of |14j . 
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renormalization group running between Z' mass scale and EW scale, which is justified 
only for a small gap between these two scale or for a low-scale Z' boson. In addition, we 
emphasize that the favored parameter regions are consistent with our assumption that 
the non-universal Z' effects in QCD penguins are negligible. This assumption requires 
lACs^sl < IAC7I or |5^J < \B^^\ < Since |5^J and \B^J\ are favored to be 

~ 10~^ and ^ 10~^ respectively, this relation can be easily accommodated. At last, to 
implement our discussions, we take a fit iii the SM and in the non-universal f/(l)' 
models for all relevant observables except ge**^. We find that the reduced value (z.e., 
X^/D.O.F.) in the SM is larger than 2, and that of the best fit in both LL and LR limits 
in the models is smaller than 1. Therefore, a better fit is obtained in the latter. 



3.2 Correlated Analysis (II) — General Case 



As discussed in Section 



and B% 



Let us 



2.1 



mixmg agam. 

Eq. (2.46) involves four of the five free parameters: \bI^^^\ and 



there are five free parameters in the general case: |-B^^ 
bcus first on 5, — 5, mixing again. The general relation in 

pj^s^- In Fig. j?! we show 
how the experimentally allowed parameter values are distributed. For each of the two 
top panels in Fig. [7} there are two peaks and two valleys toward the right. The two peaks 
in the left panel correspond to the LL limit, and the two in the right panel correspond to 
the RR limit. Meanwhile, the points in the right panel which are associated with the LL 
limit and the points in the left panel which are associated with the RR limit are localized 
in the regions ~ and |-B^| ~ 0, respectively. As for the valleys in both panels, 
they correspond to the LR limit. Observe that there are two solutions for each of these 
three limits which are specified by a difference of 180° either in or 0^ (or both). We 
only showed one of the two solutions in Fig. [T| since the difference between these two 
solutions can be resolved into B^^ as a minus sign. For the two bottom panels, the three 
special limits LL, LR and RR correspond to the bottom boundary, the diagonal line (the 
one from left-bottom to right-up) and the left boundary in the left one, respectively. 
Clearly, for the LL (RR) limit, l-B^I (|-B^|) has a relatively large value compared to the 
one in the LR limit, as seen in Fig. [TJ In the right panel, these three limits correspond 
to the two parallel bands ~ —50°, 130°, the diagonal line (the one from left-bottom 
to right-up) and the two parallel bands 0^ ~ —50°, 130°, respectively. 



In Fig. pi we illustrate how the yet free parameter B^i is constrained throug h qe 



iL,R 



using the parameter values of l-B^'"^! and (j)^^^ obtained in Fig. [?, the relation Eq. ( 2_^8 ), 
and the allowed range for qe^^ as determined in [W\ . Compare this figure with Fig. [2] we 
see that our scan selects points which are more likely to be associated with the solution 
"SI" of Bs — Bg mixing. To see the reason, let us rewrite Eq. (2.46) as 



A 



NP 



^SM 



3.59 X lO^'ilBj^fe'''''"^ + |S,^J V<; 
-2.04 X 10*^15^^,56^ |e^('^^=+'^?=) 



(3.65) 



by using the relation Eq. (1.2). Since the range of the solution "SI" is much larger 



than that of the solution "S2" at the same C.L. (see Table [5, it can be understood that 
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there are more points in the parameter space corresponding to the "SI" solution under 
the assumption of flat distribution for the random values of the relevant parameters. 

We illustrate the time-dependent CP asymmetries of the penguin-dominated modes 
in Fig. [9, taking values of \Blf'\, (pj^^^ and B^^ that are consistent with the constraints 
from Bs — Bs mixing and B^ —>■ irKs decays. As before, because of interference effects 
between the Bd — Bd mixing phase and 0^^, the points in Fig. [o] are scattered away from 
the SM limits. Hence, for each decay mode there are once again always some points 
lying in the la region. The allowed values for \B^^^ 
analysis of the Bs — Bg mixing and the Bd ^ (tt. 



and -BJ^ due to the correlated 



illustrated in Fig. [10 
the q cos ( 



, T]', p, fo)Ks CP asymmetries are 
In Fig 1]^, as done before, we map the points in Fig. 10 back to 
gsin0 plane. It is straightforward to see that the points with < —0.01 



are closer to the minima of the fit of the B 
better fit than that obtained in the SM limit. 



ttK and B ^ tth data, resulting in a 



4 Conclusions 

In this paper, we have studied the constraints on extensions of the SM with family non- 
universal U{iy gauge symmetries which result from FCNC effects in the b ^ s transi- 
tions. Using a model-independent approach in which the main requirements are family 
universal charges for the first and second generations and small fermion mixing angles, 
we have performed a correlated analysis of this set of AB = 1,2 processes. Our results 
show that within this class of models, the possible anomalies in Bg — Bg mixing and the 
time- dependent CP asymmetries of the penguin-dominated Bd — > (tt, 0, 77', p, a;, /o)-ft'5 
decays can be accommodated in a consistent way. 

Furthermore, the constraints from Bg — Bg mixing may have nontrivial implications 
not only for the hadronic decays, but also for the leptonic or semi-leptonic decays of the 



Bd mesons, as discussed in Section 2.1 As an example, recent results from BaBar 
indicate an unexpectedly large isospin asymmetry in the low dilepton mass squared 
region for combined Bd Kl^l~ and Bd — ^ K*l^l~ . In addition, K* longitudinal 
polarization and lepton forward-backward asymmetry are consistent with SM but seem 
to prefer a wrong-sign C^ operator, suggestive of right-handed currents. We leave these 
interesting issues for future exploration. 
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Appendix A Operators of Effective Hamiltonians 

A complete compilation of the relevant operators for the b ^ s transitions is given in 
the following: 

Current- Current Operators: 



Ql = (sm)v_a (m&)v-A Q2 = (SaM/3)v-A («/3^a)v-A 



(A.l) 



QCD-Penguin Operators: 



Q5 = (S&)v-A Yl 

(35 = (S&)v+A (^5) 



V-A 



V+A 



V+A 



V-A 



Q4 = (Sa6/3)v_A Y imc^Y-K 

q 

Qi = (Sa&/3)v+A Y 
Qe = (Sa6/3)v+A Y ^Woc) 



V+A 



V-A 



(A.2) 
(A.3) 
(A.4) 
(A.5) 



Electroweaik Penguin Operators: 



= ^ (S^)v-A Y ^9 (^^) 



V+A 



<37 = 2 (S&)v+A X] ^« ^^^) 

3 

= ^ (s^)v+A Y ^1 



V-A 



V+A 



= 2 (Sa&/3)v-A Y ^1 ima)y+A (^-6) 
9 

QlQ = 2 (Sa^/3)v-A X] ^9 imo^V-A (^-7) 
9 

3 

= 2 (Sa&/3)v+A Y ^9 (^/35")v-A (^-8) 
9 

3 

QlO = ^ (Sa&/3)v+A Y ^1 (^/3^a)v+A (^-9) 



Magnetic Penguin Operators: 



Q77 = ^mbSaCr'"'{l + lb)baF^,u 



QsG = ^m,s«a^^(l + i,)T-phpGl, (A.IO) 
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Semi-Leptonic Operators: 



Qgy = {bs)v-A{U)v QwA = {hs)v-A{ll) A 

Qw = (hs)v+Aill)v QwA = (hs)v+A{ll)A (A. 11) 



Bg — Bg Mixing Operators: 

Ql" = (sb)v-A{sb)v-A Q2" = {Sabi3)v-Ais/3ba)v-A 

Ql" = {sb)v+A{sb)v+A Q2" = {Sabf3)v+A{s/3ba)v+A 

Qi' = Qi' = {Sb)v+A{sb)v-A = = {sMv+A{Sf3ba)v~A (A.12) 

where indices in color singlet currents have been suppressd for simplicity, and V and A 
refer to 7^ and 7^75, respectively. 
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Appendix B Constraints of Bt{Bs 



For an order of magnitude estimate, one can temporarily ignore the effect of renor- 
malization group running. Then the branching ratio of Bg is given by (also 

see [Bl) 



Br(5, 



G 



tr — — fn mtmn 
Att 



+ 



-^\v:,Vts\' 



m 



Bs 



a 



271 sin^ 9 



w 



L uR 



BtB 



bs fi^ 



mi 



tdL tdL 



uRuR 
^^bs^^-l^J. 



{B.l] 



Here tb^ is the lifetime of Bg meson, /^^ is the corresponding decay constant, and 
Y{ml/Ml^) in the SM part is defined in [H], with 



Yix) 



X 



X 



3x 



1 



X 



l-x) 



Inx 



(B.2) 



For X = m'f/M^, we have Y{x) ~ 1. The present experimental exclusion limit at 2cr 
CL. from a combination of CDF and DO results is 



Br(5, ^ < 1.5 X 10"^ 



(B.3) 



which thus give a constraint that 
3x10-3+ ""^ 



tdL r>R 
^bs-D^,f, 



tdR tdL 
^bs-Df,f, 



+ 



nR nR 



10" 



(B.4) 



In the LR and LL limits, we have B^ 



^bs 



10 



respectively. For {Vi^ ^e^^-'^Vi^ J22 ~ 0{1) and giMz'/ig2Mz 



"3 and 



10-^ Bii = 0, 
10 — 100 (a parameter 
region favored by the correlated analysis of the anomalies in Bg —>■ Bg mixing and time- 
dependent CP asymmetries of penguin-dominated B^ meson decays; for details, see the 



last paragraph of subsection 3.1), this means that the NP contribution is comparable 
with the SM one. The experimental bound therefore can be easily satisfied. This 
conclusion also applies to the general case. 
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Appendix C Parameters 

The parameters used in our numerical analysis are summarized below: 

Mcisses, Decay Constants, Hadronic Form Factors, and Lifetimes: 



M^± = 0.139 GeV, 
Mk = 0.498 GeV, 
= 1.02 GeV, 
0.958 GeV, 
0.776 GeV, 
Mfo 
Xr, = 0.57, 
m„(/i = 4.2 GcV) = 1.86 MeV, 
m,(/i = 4.2 GcV) = 80 McV, 
rribii^ = 4.2 GcV) = 4.2 GeV, 
= 237 MeV, 
U = 130 MeV, 
= 410 McV, 
fp = 209 MeV, 
F(f^(0) = 0.330, 
Ff^{Q) = 0.379, 
Fq^{0) = 0.250, 

A^" = 0.280, 
tbo — 1.530 ps, 
Mb, = 5.37 GeV, 



M^o = 0.135 GeV, 
Mb = 5.279 GeV, 

M^ 



2.097 GcV, 
= 0.783 GcV, 
= 0.548 GeV 



0.980 GeV, 



= 0.82, 

md{^l = 4.2 GcV) = 4.22 MeV, 
mc(/i = 4.2 GcV) = 0.901 GcV, 
mtiii = Mz) = 171.7 GeV, 
fB = 190 MeV, 
fx = 160 MeV, 
= 200 McV, 
ffO_ = 180 MeV, 
""'(0) = 0.379, 
= 0.280, 
= 0.030, 

/b.V^ = 0-262 
tb- — 1.65 ps, 
TB, = 1.47 ps. 



pBKl 

%(o) 



QCD and EW Parameters: 



Gf = 1.16639 X 10"5 GeV-2, 
Mw = 80.42 GcV, 
ri2B = 0.55, 
a^iMz) = 0.118, 
A = 0.2252, 
p = 0.145, 



A% = 225 MeV, 

sin^ Ow = 0.23, 
Js = 1.627, 

a^m = 1/128, 
A = 0.8117, 
f] = 0.339, 
0.378. 



Hadronic Parameters from Lattice Calculations: 



fE 



0.262. 
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Figure 4: With i?^ and B^^ constrained by Bs — Bs mixing and B^i ttKs, the 
NP contributions to C(^^^ri' ,p,uJo)Ks ^^"^ '^{(t),rj',p,uiJo)Ks ^^'^ shown. The left 
(right) panels are the LR (LL) limits. The colors specify the C.L. that 
their inverse image points represent in Figs. [T] and [2] (yellow for la C.L. 
and blue for 2a and 1.7a C.L.). The boxes specify the experimentally 
allowed regions at Icr and 1.7a, and the dark point denotes the SM limit. 
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Figure 5: The allowed and B^^ are shown. They are constrained from 

Bs — Bs mixing (at 2a C.L.) and the fit of the B ttK (and B vrvr) 
data (at la C.L.), then selected by C(^,^',p,a;,/o)i^s' •S{<i>,rj',p,Lj,fo)Ks (at 1-7<t 
C.L. for the first four panels and 1.5a C.L. for the others). The left (right) 
panels are the LR (LL) limits. 
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Figure 6: The points in Fig. [5] are inversely mapped to the qcos(f> — qsincj) plane. 

The panels in the first row correspond to a parameter selection from 



C 



and those in the 



{<t>,r]',p,ijjo)Ksi '^{4>,v',p,uiJo)Ks ^-"^^ C.L. in Fig. |5 
second row at 1.5a C.L.. The left (right) panels are the LR (LL) limits 
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Figure 7: The distributions of \B,' \ and 



L Ft ~ 
jf^^ I a,iiLi >/)j^j resulting from Bg — Eg mixing con- 
straints are shown. The blue and purple points can be mapped to the 
experimentally allowed {Css^^f^} regions with la and 2a C.L., respec- 
tively. Here we did not distinct SI and S2 solutions any more. 
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Figure 8: The constraints on from qe^'^ are illustrated. The points of \B^^^\ 
L Fi I I 

and from Fig. 7 are randomly combined with the scattered points of 
^dd (10 '^ < \-^dd\ ^ 10~^) and then mapped to the qcos4> — qsiiK/) plane 
according to Eq. ( |2.48[ ) . The colors of the points in this plane indicate the 
C.L. that their inverse images represent in Fig. [7j 
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Figure 9: With the values of \B^^ |, and Bf^ fixed by Bg — Bs mixing and 5^ — > 

ttKs decay, the NP contributions to C(^^^rj',p,Lu,fo)Ks and Si^^^ri',p,iojo)Ks are 
illustrated in the first five panels. The colors of the points specify the 
C.L. that their inverse image points represent in Fig. [7] and Fig. |8] (yellow 
denotes la C.L. in both and blue denotes 2a and 1.7a C.L., separately). 
In the last panel, the CP asymmetries of the charmed B^ — > ipKs decay 
are presented {\Vub\ = 3.51 x 10~^ \23\)- each, the two boxes specify 
the la and 1.7a allowed regions (except for the last panel, where only the 
la box is given), and the dark point denotes the SM limit. 
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Figure 10: The distributions of </'^j'^ and are shown. The values are con- 

strained by Bs — Bg mixing {2a C.L.) and the fit of the B ^ ttK (and 
B — > tttt) data (Icr C.L.), then selected by C{^(f,^n' ,p,ll)Jo)Ks^ ^{<j>,ri' .p.u)jo)Ks 
(1.7(T C.L. for the first four panels and 1.5(T C.L. for the rest). 
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Figure 11: The points in Fig. 11 

The parameter selection from C(, 
1.7(7 C.L. in Fig. 11 the left panel, and at 1.5cj C.L. in the right panel. 



is shown at 
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